Formulae for zero-temperature conductance through a region with interaction 
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The zero-temperature linear response conductance through an interacting mesoscopic region at- 
tached to noninteracting leads is investigated. We present a set of formulae expressing the conduc- 
tance in terms of the ground-state energy or persistent currents in an auxiliary system, namely a 
ring threaded by a magnetic flux and containing the correlated electron region. We flrst derive the 
conductance formulae for the noninteracting case and then give arguments why the formalism is 
also correct in the interacting case if the ground state of a system exhibits Fermi liquid properties. 
We prove that in such systems, the ground-state energy is a universal function of the magnetic flux, 
where the conductance is the only parameter. The method is tested by comparing its predictions 
with exact results and results of other methods for problems such as the transport through sin- 
gle and double quantum dots containing interacting electrons. The comparisons show an excellent 
quantitative agreement. 



I. INTRODUCTION 

The measurements of the conductivity and the elec- 
tron transport in general are one of the most direct and 
sensitive probes in solid state physics. In such measure- 
ments many interesting new phenomena were signaled, 
in particular superconductivity, transport in metals with 
embedded magnetic impurities and the related Kondo 
physics, heavy fermion phenomena and the physics of the 
Mott-Hubbard transition regime. In the last decade tech- 
nological advances enabled controlled fabrication of small 
regions connected to leads and the conductance, relating 
the current through such a region to the voltage applied 
between the leads, also proved to be a relevant property 
of such systems. There is a number of such examples, e.g. 
metallic islands prepared by e-beam lithography or small 
metallic grains^, semiconductor quantum dots^, or a sin- 
gle large molecule such as a carbon nanotube or DNA. It 
is possible to break a metallic contact and measure the 
transport properties of an atomic-size bridge that forms 
in the break^ , or even measure the conductance of a single 
hydrogen molecule, as reported recently in Ref.". In all 
such systems, strong electron correlations are expected 
to play an important role. 

The transport in noninteracting mesoscopic systems 
is theoretically well described in the framework of the 
Landauer-Biittiker formalism. The conductance is deter- 
mined with the Landauer-Biittiker formula^"^, where the 
key quantity is the single particle transmission amplitude 
t(£) for electrons in the vicinity of the Fermi energy. The 
formula proved to be very useful and reliable, as long as 
electron-electron interaction in a sample is negligible. 

Although the Landauer-Biittiker formalism provides a 
general description of the electron transport in noninter- 
acting systems, it normally cannot be used if the interac- 
tion between electrons plays an important role. Several 
approaches have been developed to allow one to treat also 
such systems. First of all, the Kubo formalism provides 
us with a conductance formula which is applicable in the 
linear response regime and has, for example, been used 



to calculate the conductance in Refs.*'®. A much more 
general approach was developed by Meir and Wingreen 
in Ref.^°. Within the Keldysh formalism they manage 
to express the conductance in terms of nonequilibrium 
Green's functions for the sample part of the system. The 
formalism can be used to treat systems at a finite source- 
drain voltage and can also be extended to describe time- 
dependent transport phenomena^ ^. The main theoretical 
challenge in these approaches is to calculate the Green's 
function of a system. Except in some rare cases where 
exact results are available, perturbative approaches or 
numerical renormalization group studies are employed. 

In this paper we propose an alternative method for 
calculating the conductance through such correlated sys- 
tems. The method is appHcable only to a certain class of 
systems, namely to those exhibiting Fermi liquid prop- 
erties, at zero temperature and in the linear response 
regime. However, in this quite restrictive domain of va- 
lidity, the method promises to be easier to use than the 
methods mentioned above. We show that the ground- 
state energy of an auxiliary system, formed by connecting 
the leads of the original system into a ring and threaded 
by a magnetic flux, provides us with enough informa- 
tion to determine the conductance. The main advantage 
of this method is the fact that it is often much easier 
to calculate the ground-state energy (for example, using 
variational methods) than the Green's function, which 
is needed in the Kubo and Keldysh approaches. The 
conductance of a Hubbard chain connected to leads was 
studied recently using a special case of our method and 
DMRG^^'^^ and a special case of our approach was ap- 
plied in the Hartree-Fock analysis of anomalies in the 
conductance of quantum point contacts^". The method 
is related to the study of the charge stiffness and persis- 
tent currents in one-dimensional systems^^"^*. 

The paper is organized as follows. In Section II we 
present the model Hamiltonian for which the method is 
applicable. In Section III we derive general formulae for 
the zero-temperature conductance through a mesoscopic 
region with noninteracting electrons connected to leads. 
In Section IV we extend the formalism to the case of in- 
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teracting electrons. We give arguments why the formal- 
ism is correct as long as the ground state of the system 
exhibits Fermi hquid properties. In Section V conver- 
gence tests for a typical noninteracting system are first 
presented. Then we support our formalism also with nu- 
merical results for the conductance of some non-trivial 
problems, such as the transport through single and dou- 
ble quantum dots containing interacting electrons and 
connected to noninteracting leads. These comparisons, 
including the comparison with the exact results for the 
Anderson model, demonstrate a good quantitative agree- 
ment. After the conclusions in Section VI we present 
some more technical details in Appendix A. In Appendix 
B we describe the numerical method used in Section V. 



II. MODEL HAMILTONIAN 
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Figure 1: Schematic picture of the system described by Hamil- 
tonian (1). 

In this Section we introduce a general Hamiltonian de- 
scribing a mesoscopic sample coupled to leads as shown 
in Fig. 1. The Hamiltonian is a generalization of the well 
known Anderson impurity model^^. We split the Hamil- 
tonian into five pieces 

H = HL+VL + Hc + VTi + HR, (1) 

where He models the central region, Hl and Hr describe 
the left and the right lead, and Vl and Vr are the tunnel- 
ing couplings between the leads and the central region. 
We can also split the Hamiltonian into a term H^^^ de- 
scribing independent electrons and a term U describing 
the Coulomb interaction between them 

H = + U. (2) 

One can often neglect the interaction in the leads and 
between the sample and the leads. We assume this is 
the case. Then the central region is the only part of the 
system where one must take the interaction into account 

He = H^c^ + U. (3) 
Here H^^^ describes a set of noninteracting levels 

4°^ = E H^Ad^^' (4) 



where d]^ (dia) creates (destroys) an electron with spin 
a in the state i. The states introduced here can have 
various physical meanings. They could represent the 
true single-electron states of the sample, for example dif- 
ferent energy levels of a multi-level quantum dot or a 
molecule. In this case, the matrix H^^^ is diagonal and 
its elements are the single-electron energies of the sys- 
tem. Another possible interpretation of Hamiltonian (4) 
is that the states i are local orbitals at different sites of 
the system. In this case, the diagonal matrix elements 
of i?,?'' are the on-site energies for these sites, while the 
off-diagonal matrix elements describe the coupling be- 
tween different sites of the system. The sites could have 
a direct physical interpretation, such as dots in a dou- 
ble quantum dot system or atoms in a molecule, or they 
could represent fictitious sites obtained by discretization 
of a continuous system. There are other possible choices 
of basis states for the central region. For example, in a 
system consisting of two multi-level quantum dots one 
could use single-electron basis states for each of the dots 
and describe the coupling between the dots with tunnel- 
ing matrix elements. 

The Coulomb interaction between electrons in the sam- 
ple is given by an extended Hubbard-type couphng 

where Uic = d\^dic is the operator counting the number 
of electrons with spin a at site i. For convenience, we 
wrote down only the expression for the Coulomb inter- 
action in the case, where basis states represent different 
sites in real space. The expression becomes more com- 
plicated if a more general basis set is used. 

We describe the leads or contacts as two semi-infinite, 
tight-binding chains 

Hi^(R)=-to E c\^Ci+ia+^-c., (6) 

i,i+leI-(H) 

where cj^ (cjcr) creates (destroys) an electron with spin 
a on site i and to is the hopping matrix element between 
neighboring sites. Such a model adequately, at least for 
energies low or comparable to to, describes a noninteract- 
ing, single-mode and homogeneous lead. It would be easy 
to generalize the lead Hamiltonian to describe a more 
realistic system, for example by modehng the true ge- 
ometry or ahowing for a self-consistent potential due to 
interaction between electrons. However, the physics we 
are interested in, is usually not changed dramatically by 
not including these details into the model Hamiltonian 
and therefore, we will not discuss this issue into detail. 

Finally, there is a term describing the coupling between 
the sample and the leads, 

VL(fi) = E yL(R)jic]„dia + h..C., (7) 
iec 
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where VL(/?)ij is the hopping matrix element between 
state i in the sample and site j in a lead. 
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Figure 2: The sample embedded in a ring formed by joining 
the left and right leads of the system in Fig. 1. Magnetic flux 
3> penetrates the ring. 

In the following Sections wc discuss the conductance 
through the system introduced above. To derive the con- 
ductance formulae, we will need a slightly modified sys- 
tem. This auxiliary system is a ring formed by connecting 
the ends of the left and right leads of the original system 
as shown in Fig. 2. The ring is threaded by a magnetic 
flux <I> in such a way that there is no magnetic field in the 
region where electrons move. We can then perform the 
standard Peierls substitution^^ and transform the hop- 
ping matrix elements of the Hamiltonian (1) according 
to 

i..^i..eif^7A.dx^ (8) 

where Xj is the position of site i and A is the vector 
potential due to the flux, obeying 



$ = -0 = ^ A • dx. (9) 

Here wc defined a dimensionlcss magnetic fiux (t>. The 
energy of the system is periodic in 4> with a period of 27r 
and depends only on value of (/) and not on any details of 
how the flux is produced. If the original Hamiltonian (1) 
obeys the time-reversal symmetry, the energy does not 
change if the magnetic field is reversed, 

E {-<!>) =E{cl>). (10) 

III. CONDUCTANCE OF A 
NONINTERACTING SYSTEM 

In this Section we limit the discussion to noninteract- 
ing systems, i.e. we set J7 = in Eq. (2). In such systems, 
the Landauer-Biittiker formula^"^ 

G = Go\t{eF)\\ (11) 



which relates the zero-temperature conductance G to 
the transmission probability \t (£f)|^ for electrons at the 
Fermi energy ep, can be appHed. The proportionality 
coefficient, Gq = is the quantum of conductance. 
Below we first derive a set of formulae, which relate the 
transmission probability, and consequently the conduc- 
tance, to single-electron energy levels of the auxiliary ring 
system introduced in the previous Section. Then we de- 
rive another set of formulae, relating the conductance to 
the ground-state energy of the auxifiary system. One of 
these formulae was derived before in Ref.^'', and a limit- 
ing case of another one was discussed in Refs.^^'^^. Here 
we present a unified approach to the problem, from which 
these results emerge as special cases. 

A. Formulae relating conductance to 
single-electron energy levels 

Let us consider eigenstates of an electron moving on a 
ring system introduced in the previous Section. We will 
be interested only in energies of these states and not in 
the precise form of wavcfunctions. The energy of an elec- 
tron on a ring penetrated by a magnetic fiux (j) depends 
only on the magnitude of the fiux and therefore, any vec- 
tor potential fulfilling condition (9) is good for our pur- 
pose. We choose a vector potential constant everywhere 
except between sites io and io + l, both in the lead part of 
the ring as shown in Fig. 2. The hopping matrix element 
between the two sites is thus modified to toe^'^. With 
no fiux penetrating the ring, the electron's wavefunction 
in the lead part of the system is ae**^* -|- 6e~*'^*, where k 
is the electron's wavevector and a and b are amplitudes 
determined by properties of the central region. If there 
is a fiux through the ring, the wavefunction is modified. 
The Schrodinger equations for sites io and io + l show 
us that the appropriate form is ae**^* + he^^^"^ for i < io 
and ae~*'^e*'°* -|- 6e~*'^e~*'^* for i > io. The scattering 
matrix of the central region provides a relation between 
coefficients a and b, 

( 6e-*'^e^'=^ \ _( rk f'k\( ae-'^e-'>'^ \ 

a )-[tk r'j[ b )■ (12) 

The elements of the scattering matrix, tk and rk {t'/, and 
r'l^), arc the transmission and reflection amplitudes for 
electrons coming from the left (right) lead, and N is the 
number of sites in the lead part of the ring. We added 
phase factors e±*'=^ to the "left lead" amplitudes to com- 
pensate for the phase difference an electron accumulates 
as it travels through the lead part of the ring. The scat- 
tering matrix defined this way does not depend on N 
and (j), and equals the scattering matrix of the original, 
two-lead system. Eq. (12) is a homogeneous system of 
linear equations, solvable only if the determinant is zero. 
Using the unitarity property of the scattering matrix, the 
eigenenergy condition becomes 

i'fee'^ + tke-''^ = e''^^ + %e-''=^. (13) 
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We assume the Hamiltonian of the original, two lead sys- 
tem obeys the time-reversal symmetry and therefore, the 
scattering matrix is symmetric^^, tk = t'^. Expressing 
the transmission amplitude in terms of its absolute value 
and the scattering phase shift tk — \tk \ e"^'', we arrive at 
the final form of the eigenenergy equation 

\t k\ cos (p ^ cos [kN - 'fk) ■ (14) 

In Fig. 3 a graphical representation of this equation is 
presented. 




Figure 3: A graphical representation of the eigenvalue equa- 
tion (14). The shaded region represents the allowed values 
of the left hand side of the equation for different values of 
magnetic flux (for example, the dashed line shows the values 
for (j) — j). The full line represents the right hand side of 
the eigenvalue equation. The system is presented in Fig. 7, 
N = 100. 

To extract the transmission probability \tkf , we pro- 
ceed by differentiating the eigenvalue equation with re- 
spect to cos (p 

d\tk\ 



dcos( 



■cos^-l- \tk\ 



= - sin {kN + ipk)(N 



dk 



dcos ( 



d(pk 
d cos ( 



= ±Vl-|*.rcos^4A^^ + /^). (15) 
* \ <J cos (j) o cos (j) J 

The sign of the last expression depends on weather k 
belongs to a decreasing (+) or an increasing (— ) branch 
of the cosine function in Eq. (14), or equivalently, if we 
are interested in an eigenstate with odd (-I-) or even (— ) 
n, where n indexes the cigcnstates from the one with the 
lowest energy upward. Let us choose an eigenstate and 
consider how the corresponding wavevector k changes as 
the magnetic flux (f) is varied from to tt. It is evident 
that the variation in k is of the order of as the cosine 
function in the right hand side of Eq. (14) oscillates with 
such a period. Let as assume that the number of sites 
in the ring is large enough that transmission amplitude 
does not change appreciably in this interval 



dtk 
dk 



(16) 



Then the derivatives 



dk 



d\tk\ 



and 



order of and Eq. (15) simplifies to 
\tk\ = ±^l-|ife|^cos2 ^N- 



dcos 4 



are of the 



dk 



d cos 



O 



(17) 



Introducing the density of states in the leads p (e) = 
which for example, for a tight-binding lead with 
only nearest-neighbor h opping and dispersion £k = 
—2tocosk equals l/(7rA/4t^^^e|), we finally obtain 



d arccos (=F | i (^fc ) | cos < 
dcoscj) 



= irNp (£fe) 



dsk 
9 cose 



(18) 



where t {sk) = tk- The condition Eq. (16) of validity can 
also be expressed in a form involving energy as a variable 



1 



Pie) 



dt (e) 



de 



(19) 



Eq. (18) is the central result of this work. It expresses the 
transmission probability |i(£)|^ of a sample connected 
to two leads in terms of the variation of single-electron 
energy levels with magnetic flux penetrating the auxiliary 
ring system. Employing the Landauer-Biittiker formula 
Eq. (11), this result also provides the zero-temperature 
conductance of the system. From the derivation it is 
evident that the method becomes exact as we approach 
the thermodynamic limit N ^ oo. 

In general Eq. (18) has to be solved numerically to 
obtain the transmission probability on a discrete set of 
energy points, one for each energy level of a system. By 
increasing the system size N, the density of these points 
increases and the errors decrease, as the condition (19) is 
fulfllled better. We will return to this point in Section V 
where we consider the convergence issues in detail. Here 
we present some special cases of Eq. (18) where analytic 
expressions can be obtained. By averaging the equation 
over values of flux (j) between ^ = and (j) = n (note 
that we may treat |i(£fe)| and p{£k) as constant while 
averaging as the resulting error is of the order of ;^), 
we can relate the transmission probability to the average 
magnitude of the derivative of a single-electron energy 
with respect to the flux: 



\t{ek) 



sm 



-Np [Ek) 



dsk 



(20) 



Note that it is enough to calculate the energy levels at 
= a nd (p = T T to calculate the transmission proba- 
bility as \dek/d^\ = \ \ek (tt) - Sk (0)|. In Fig. 4(a) it is 
illustrated how a large variation of single-electron energy 
as the flux is changed from (/) = to (i> = tt corresponds 
to a large conductance and vice versa. The transmission 
probability can also be calculated from the derivative at 
<^='^ resulting in the second formula 



\i(ek)\'= UArp(£fc) 



dsk 



(21) 
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Again, Fig. 4(a) shows that there is a correspondence 
between a large sensitivity of a single-electron energy to 
the flux at 4> = ^, and a large conductance. Finally, we 
observe that the curvature of energy levels aX (j) = and 
(p = IT also gives information of conductance. The third 
formula reads 



(22) 



sum of single-electron energies up to the Fermi energy 
ep, multiphed by 2 because of the electron's spin 



E = 2 



(23) 



We will show that the transmission probability at the 
Fermi energy |t(eF)| is related to the ground-state en- 
ergy of the ring system 




Figure 4: (a) Single-electron energy levels (full lines) and the 
ground-state energies when a given single-electron level is at 
the Fermi energy (dashed lines). Ground-state energies are 
shifted so that both curves coincide for (j) = ^. Note that 
the energy curves are symmetric about (f> = as required by 
Eq. (10). The shaded area represents the magnitude of the 
transmission amplitude. The system and the energy interval 
is the same as in Fig. 3. (b, c) The large N universal form of 
the ground-state energy vs. flux curve for an even (b) and an 
odd (c) number of electrons in a system. The magnitude of 
the transmission amplitude goes from to 1 in steps of 0.1. 



B. Formulae relating conductance to the 
ground-state energy 

Above we showed how the flux variation of the en- 
ergy of the last occupied single-electron state allows one 
to calculate the zero-temperature conductance through 
a nonintcracting sample. The goal of this section is to 
derive an alternative set of formulae, expressing the zero- 
temperature conductance in terms of the flux variation 
of the ground-state energy E, which for an even num- 
ber of electrons in a noninteracting system is simply a 



1 Sarccos^ (=F \t (£f) \ cose 

TT 



9cos( 



ttNp (sf) 



dE 

d cos (f) ' 



(24) 



where the sign is — and -f- for an odd and an even num- 
ber of occupied single-electron states, respectively. The 
expression Eq. (24) is evidently correct if there are no 
electrons in the system, as it gives a zero conductance in 
this case. To prove the formula for other values of the 
Fermi energy, we use the principle of the mathematical 
induction. To simplify the notation we introduce 



fs (en) = N 



d£n 

9 cose 



1 d arccos {s\t (£„ ) | cos ( 



(En) 



9cos< 



(25) 



Fs (£n) = N 



dE„ 



1 d arccos^ {s\t (£„) | cos 0) 



dcosi 



(26) 

where En is the ground-state energy of a system with 
the Fermi energy at £„ and s is either 1 or —1, depend- 
ing on the signs in Eqs. (18) and (24). Differentiating the 
relation En = En-i + '2.£n with respect to cos0, express- 
ing the result in terms of functions fs and Eg introduced 
above, and making use of the fact that the sign s alter- 
nates with n, we obtain 



Fs{en) = F-s (£n-l) + 2/,(£„). 



(27) 



If we manage to show that this really is an identity, we 
have a proof of Eq. (24) . Using the exact relation Fs (e) — 
F-s (e) = 2/s (e), the expression transforms into 

Fs (e„) = Fs (e„) - [F.s (£„) - F-,, (£„-i)] • (28) 

For a large number of sites N in the ring and corre- 
spondingly, a small separation of single-electron energy 
levels which is of the order of , the term in parenthe- 
sis equals F'_g (£„) (e„ — e„_i). F'_g (£„) can be factored 
into sF (e„) where F (e„) does not depend on sign s. 
Therefore, although the term in parenthesis in of the or- 
der of its sign alternates for successive energy levels 
while its amplitude stays the same. Thus the error in- 
duced by this term does not accumulate, it just adds an 
additional error of the order of to the flnal result. 

In Fig. 4(a), the variation of the ground-state energies 
with magnetic flux is compared to the variation of the 
corresponding single-electron energies. Note that as a 
consequence of Eq. (24), the ground-state energy in the 
large N limit takes a universal form (see Fig. 4(b)) 
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E{cI>)-eQ = 



( arccos^ \t {e f) \ cos c 



For systems with an odd number of electrons, the ground- 
state energy is obtained by adding a single-electron en- 
ergy corresponding to Eq. (18) and the universal form 
reads (see Fig. 4(c)) 



E{<l>)-E (I) = ^2^1 (sf) ^'''^'''^ (£f)| cos<^) . 

(30) 

In general, Eq. (24) can only be solved numerically 
to obtain the transmission probability. However, as was 
the case for single-electron energies, analytic solutions 
can be found in certain special cases. The derivative of 
the ground-state energy with respect to flux gives the 
persistent current in the ring j = f §f Using the 

Landauer-Biittiker formula Eq. (11), one can calculate 
the conductance from the flux averaged magnitude of the 
persistent current in the system 



\t{eFr=sm'(—Np{eF)\jm 



(31) 



Only two ground-state energy calculations n eed to be 
performed to obtain the conductance as f|.j((/')| = 
^ \E (tt) — E {0)\. This formula was also discussed in 
Refs.^^'^^ for the case where the transmission probability 
is small. The second formula relates the conductance to 
the persistent current at 



TT 14,24 
2 ' 



tieF)f=r-^NpieF)j{l) 



(32) 



The third formula, corresponding to Eq. (22) in the 
single-electron case, turns out to be more complicated 
and gives an implicit relation for \t (£f)| 



ttNp (sf) 



±- 



d^E 



^ mm, max 

2\t{eF)\ 



arccos(± |t (eir)l) . (33) 



nJl-ltisFW 



Here the upper and the lower signs correspond to the sec- 
ond derivative at a minimum and at a maximum of the 
energy vs. flux curve, respectively. Minima (maxima) 
occur at ^ = (tt) if an odd number of single-electron 
levels is occupied and at = tt (0) if an even number of 
levels is occupied. The second derivative in a minimum is 
proportional to the charge stiffness D — ^ '^^^/^^^Imin 
of the system^^'^^. We can also define the corresponding 
quantity for a maximum as D = — d'^E/d(P\ . In 

^ ^ z / ■ I max 

general, Eq. (33) has to be solved numerically. However, 



in the limit of a very small conductance and in the vicin- 
ity of the unitary limit, additional analytic formulae are 
valid 



'*^'^^'-U + ¥27rp(.F)A 



\t{eF)\ 
\t{eF)\ 



•0, 
. 1. 



(34) 



Note that there is a quadratic relation between the con- 
ductance and the charge stiffness in the low conductance 
limit. The corresponding formulae for the maximum of 
the energy vs. fiux curve are 



\t{SF) 



27rp {sf)D, 



(2-kp(sf)D) 



\t{eF)\ 

\t{eF)\ 



0, 

1. 



(35) 



A detailed analysis of convergence properties of the for- 
mulae derived in this Section is presented Section V. 

We stress again that the validity of these formulae 
is based on an assumption that the number of sites in 
the ring is sufficiently large according to the condition 
Eq. (19). This means that if t{e) exhibits sharp reso- 
nances, the calculation has to be performed on such a 
large auxiliary ring system that in the energy interval 
of interest (the width of the resonance) there is a large 
number of eigenenergies £„. Then t{£) <^ t{en'), where 
£„' is the eigenenergy closest to e. Such sharp resonances 
in t{e) are expected e.g. in chaotic systems^^'^''. The 
present method might be impractical (but still correct) 
in this case. 



IV. CONDUCTANCE OF AN INTERACTING 
SYSTEM 

The zero-temperature conductance of a noninteract- 
ing system can thus be determined with the transmission 
probability obtained from one of the formulae we derived 
in the previous Section, and the Landauer-Biittiker for- 
mula. The main challenge, however, remains the ques- 
tion of the validity of this type of approach for interact- 
ing systems. In this Section we give arguments why the 
approach is correct for a class of interacting systems ex- 
hibiting Fermi liquid properties. In order to reach this 
goal, we present four essential steps as follows. 



Step 1: Conductance of a Fermi liquid system at 
T — 



The basic property that characterizes Fermi liquid 
systcms^^ is that the states of a noninteracting system of 
electrons are continuously transformed into states of the 
interacting system as the interaction strength increases 
from zero to its actual value. One can then study the 
properties of such a system by means of the perturba- 
tion theory, regarding the interaction strength as the per- 
turbation parameter. The concept of the Fermi liquid 
was first introduced for translation-invariant systems by 
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Landau^^'^°, and was later also extended to systems of 
the type we study here^^. 

The linear response conductance of a general interact- 
ing system of the type shown in Fig. 1 can be calculated 
from the Kubo formula^'^^ 



le 



G = lim -11// (u + iS) , 

w^O LO + 10 



(36) 



where 11// {u + iS) is the retarded current-current corre- 
lation function 

H// {t - t') = -i9 {t - t') {[I {t) , I {t')]) . (37) 

For Fermi liquid systems at T = 0, the current-current 

correlation function can be expressed in terms of the 
Green's function G„'„ (z) of the system and the conduc- 
tance is given with^ 



2e2 



1 



-inp (ef) 



(38) 



where n and n' are sites in the left and the right lead, 
respectively. One can define the transmission amplitude 
as 



t{s) 



1 



-inp (e 



■e-''=("'-")G„.„(e + z5), 



and the conductance formula Eq. (38) then reads 



(39) 



(40) 



For non-interacting systems, t (e) defined this way re- 
duces to the standard transmission amplitude (Fisher- 
Lee relation^^) and Eq. (40) represents the Landauer- 
Biittiker formula. In the next step, we will show that 
the transmission amplitude Eq. (39) has a direct physi- 
cal interpretation also for interacting systems, being the 
transmission amplitude of Fermi liquid quasiparticles. 



Step 2: Quasiparticle Hamiltonian 

In this step, we generalize the quasiparticle approxi- 
mation to the Green's function, presented for the single- 
impurity Anderson model in Ref.^'', to the case where the 
interaction is present in more than a single site. 

In Fermi liquid systems obeying the time-reversal sym- 
metry, the imaginary part of the retarded self-energy at 
T = vanishes at the Fermi energy and is quadratic 
for frequencies close to the Fermi energy^^'^^. Using 
the Fermi energy as the origin of the energy scale, i.e. 
(J — ef — * w, we can express this as 



ImS {lo + iS) oc Lo'^ 



(41) 



Close to the Fermi energy, the self-energy can be ex- 
panded in powers of lu resulting in an approximation to 



the Green's function, 

G"^ {uj + id) =oj1- H(°) --EiO + id)- (42) 
+ O . (43) 



—OJ 



as (w + i5) 



duj 



Here H*^"^ contains matrix elements of the noninteracting 
part of the Hamiltonian (2). Note that expansion coeffi- 
cients are real because of Eq. (41). Let us introduce the 
renormalization factor matrix Z as 



z-i = i- 



5S {io + id) 



(44) 



w=0 



The Green's function for u) close to the Fermi energy can 
then be expressed as 

(w + iS) = (w + iS) Z-V2 + (oj^) , (45) 

where we defined the quasiparticle Green's function 

G"^ (w + iS)=ojl-il (46) 

as the Green's function of a noninteracting quasiparticle 
Hamiltonian 



H = ZV2 



H(°) -h S (0 -h iS) 



Zi/2. 



(47) 



Note that factoring the renormalization factor matrix 
as we did above ensures the hermiticity of the resulting 
quasiparticle Hamiltonian. 

Matrix elements of Z differ from those of an identity 
matrix only if they correspond to sites of the central re- 
gion. In other cases, as the interaction is limited to the 
central region, the corresponding self-energy matrix el- 
ement is zero. Therefore, comparing the quasiparticle 
Hamiltonian to the noninteracting part of the real Hamil- 
tonian, we observe that the effect of the interaction is 
to renormalize the matrix elements of the central region 
Hamiltonian (4) and those corresponding to the hopping 
between the central region and the leads (7). The values 
of the renormalized matrix elements depend on the value 
of the Fermi energy of the system. 

Let us illustrate the ideas introduced above for the 
case of the standard Anderson impurity model'^'' . We cal- 
culated the self-energy in the second-order perturbation 
theory approximation^'''"^^ and constructed the quasipar- 
ticle Hamiltonian according to Eq. (47). In Fig. 5 the lo- 
cal spectral functions corresponding to both the original 
interacting Hamiltonian and the noninteracting quasipar- 
ticle Hamiltonian are presented. The agreement of both 
results is perfect in the vicinity of the Fermi energy where 
the expansion (42) is valid. 

The reason for introducing the quasiparticle Hamilto- 
nian is to obtain an alternative expression for the con- 
ductance in terms of the quasiparticle Green's function. 
Eq. (45) relates the values of the true and the quasipar- 
ticle Green's function at the Fermi energy. 



G(0 + i(5) = Zi/2G(0 + i(5)Zi/2. 



(48) 
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Figure 5: The T = local spectral function and the corre- 
sponding quasiparticle approximation for the Anderson im- 
purity model as shown in Fig. 9. The values of parameters 
are ti — OAto, U = 1.92io and Sd = — -j. The calculations 
were performed within the second-order perturbation theory 
as described in Appendix A. 



are sites in the leads, 
as a consequence of the 



Specifically, if both n and n 

Gn'n{0 + i6) = Gn'n{0 + tS) 

properties of the renormalization factor matrix Z dis- 
cussed above. Eq. (39) then tells us that the zero- 
temperature conductance of a Fermi liquid system is 
identical to the zero-temperature conductance of a non- 
interacting system defined with the quasiparticle Hamil- 
tonian for a given value of the Fermi energy. 



Step 3: Quasiparticles in a finite system 

The conclusions reached in the first two steps are based 
on an assumption of the thermodynamic limit, i.e. tlicy 
are valid if the central region is coupled to semiinfinite 
leads. Here we generahze the concept of quasiparticles 
to a finite ring system with N sites and M electrons, 
threaded by a magnetic fiux 0. Let us define the quasi- 
particle Hamiltonian for such a system, 

H (iV, 0; M) = rjj(o) 0) _^ s (0 id)] Z^/^ 

(49) 

Here the self-energy and the renormalization factor ma- 
trix are determined in the thermodynamic limit where, 
as we prove in Appendix A, they are independent of (j) 
and correspond to those of an infinite two-lead system. 

Suppose now that we knew the exact values of the 
renormalized matrix elements in the quasiparticle Hamil- 
tonian (49). As this is a noninteracting Hamiltonian, we 
could then apply the conductance formulae of the pre- 
vious Section to calculate the zero-temperature conduc- 
tance of an infinite two-lead system with the same central 
region and central region- lead hopping matrix elements, 
i.e. of a system described with the quasiparticle Hamil- 
tonian (47). As shown in step 2, this procedure would 



provide us with the exact conductance of the original in- 
teracting system. However, to obtain the values of the 
renormalized matrix elements, one needs to calculate the 
self-energy of the system, which is a difficult many-body 
problem. In the next step, we will show, that there is an 
alternative and easier way to achieve the same goal. 

In Fig. 6 we compare the spectral density of an Ander- 
son impurity embedded in a finite ring system to that of 
the corresponding quasiparticle Hamiltonian (49). Note 
that the spectral density of the quasiparticle Hamilto- 
nian correctly describes the true spectral density in the 
vicinity of the Fermi energy. 




-0.15 



0.15 



Figure 6: (a) The T = local spectral function as in Fig. 5, 
but for a ring system with A'" = 400 sites and flux = 
(b) The spectral function in the vicinity of the Fermi energy 
(dashed lines) compared to that corresponding to the quasi- 
particle Hamiltonian (49). Both the spectral density of the 
interacting system and the matrix elements of the quasiparti- 
cle Hamiltonian were calculated within the second order per- 
turbation theory. 



Step 4: Validity of the conductance formulae 

In this last step we finally show how to calculate the 
conductance of an interacting system. In Appendix A 
we study the excitation spectrum of a finite ring system 
threaded with a magnetic fiux and containing a region 
with interaction. We show that 



E[N,(I);M + 1]- E [N, 0; M] = 
= e{N,cl>-M-l) + o[N-i^ 



(50) 



where E [N. (f): M) and E {N, &, M + 1) are the ground- 
state energies of the interacting Hamiltonian for a ring 
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system with N sites and flux 4>, containing M and M + 1 
electrons, respectively, and e (TV, 0; M; 1) is the energy 
of the first single-electron level above the Fermi en- 
ergy of the finite ring quasiparticle Hamiltonian (49). 
This estimation allows one to use single-electron for- 
mulae of Sec. Ill A to calculate the zero-temperature 
conductance for a Fermi liquid system. We showed in 
step 3 that inserting s {N, </>; M; 1) into these formu- 
lae would give us the correct conductance. Eq. (50) 
proves, that the same result is obtained if the differ- 
ence of the ground state energies of an interacting system 
E [N, (l);M +1]-E [TV, (f): M] is inserted into a formula in- 
stead. The estimated error, which is of the order of 2 , 
is for a large N negligible, because it is much smaller than 
the quasiparticle level spacing, which is of the order of 

As demonstrated in Sec. IIIB, the conductance of a 
noninteracting system can also be calculated from the 
variation of the ground-state energy with fiux through 
the ring. The proof of the formulae involved only the 
properties of a set of neighboring single-electron energy 
levels. We assumed the validity of single-electron con- 
ductance formulae for each of these levels and made use 
of the fact that the ground-state energy of the system 
increases by a sum of the relevant single-electron ener- 
gies as the levels become occupied with additional elec- 
trons. For Fermi Hquid systems, the first assumption 
was proved above. The second assumption, which for 
noninteracting sj^stems is obvious, is proved in Appendix 
A. There we show that as a finite number of additional 
electrons is added to an interacting system, the succes- 
sive ground-state energies are determined by the single- 
electron energy levels of the same quasiparticle Hamil- 
tonian with a very good accuracy (Al). Therefore, the 
proof of Sec. IIIB is also valid for interacting Fermi Hq- 
uid systems, provided a system is the Fermi liquid for all 
values of the Fermi energy below its actual value. 

V. NUMERICAL TESTS OF THE METHOD 
A. Noninteracting system 




Figure 7: A double barrier noninteracting system. The height 
of the barriers is 0.5to, where to is the hopping matrix element 
between neighboring sites. 

In this Section wc discuss the convergence properties 
of the conductance formulae derived in Sec. III. As a 
test system we use a double-barrier potential scattering 
problem presented in Fig. 7. Results of various formulae 
for different number of sites in the ring are presented in 



Fig. 8. The exact zero-temperature conductance for this 
system exhibits a sharp resonance peak superimposed on 
a smooth background conductance. We notice immedi- 
ately that as the number of sites in the ring increases, 
the convergence is generally faster in the region where 
the conductance is smooth than in the resonance region, 
which is consistent with the condition (19). Comparing 
the results obtained employing different conductance for- 
mulae we observe that the convergence is the fastest in 
both the single-electron and the ground-state energy case 
if the formulae of Eqs. (22) and (33) are appHed to the 
maximum of the energy vs. flux curve (or to the mini- 
mum in the single-electron case). Formulae of Eqs. (20) 
and (31) expressing the conductance in terms of the dif- 
ference of the energies at = and = tt converge 
somewhat slower. Note however that in the former case 
the second derivative of the energy with respect to the 
flux has to be evaluated while in the later, the energy 
difference is large and because of that, the calculation is 
much more well behaved. From the computational point 
of view there is another advantage of the energy differ- 
ence formulae. In this case, all the matrix elements can 
be made real if one chooses such a vector potential that 
only one hopping matrix clement if modified by the flux 
as then the additional phase factor is e^^"^ = — 1. Finally, 
the remaining formulae, employing the slope of the en- 
ergy vs. flux curve at <^ = J and the curvature in the 
minimum of the ground-state energy vs. flux curve, do 
not show convergence properties comparable to those of 
the formulae discussed above. 



B. Anderson impurity model 

In 1980s several theories^°'^^ were put forward propos- 
ing a realization of the Anderson impurity model^*' in 
systems consisting of a quantum dot coupled to two leads 
(see Fig. 9). These theories show that the topmost occu- 
pied energy level in a quantum dot with an odd number 
of electrons can be associated with the Anderson model 
£d level and such a system should mimic the old Kondo 
problem of a magnetic spin | impurity in a metal host. In 
recent years signatures of the Kondo physics in electron 
transport through quantum dots have also been found 
experimentally^^'^^ . The Anderson model is well defined 
and is an attractive testing ground for new numerical 
and analytical methods that are developed to tackle other 
challenging many-body problems. Therefore, we will also 
take it as a nontrivial example to test results of the con- 
ductance formulae we derived in this paper. 

There are three distinct parameter regimes of the An- 
derson model. If Ed < £f < £d + U with \ed + U — £f| ^ 
A and \ed — sf\ ^ A, where A is the coupling of the 
quantum dot to leads, wc arc in the Kondo regime. In this 
regime, a narrow Kondo resonance is formed in the spec- 
tral function at the Fermi energy for temperatures below 
and close to the Kondo temperature, which corresponds 
to the width of the resonance. The zero-temperature 
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Figure 8: Exact and approximate zero-temperature conductance vs. Fermi energy curves for the system in Fig. 7. The shaded 
area shows the exact result. The left set of figures shows the approximations obtained using the single-electron formulae of 
Sec. Ill A while the right set of figures corresponds to the ground-state energy formulae of Sec. IIIB. Different curves correspond 
to different number of sites in the ring. In figures (ai) and (a2) the conductance was calculated using Eqs. (20) and (31), in 
figures (bi) and (1)2) using Eqs. (21) and (32), while in the other figures Eqs. (22) and (33) were used, in (ci) and (02) applied 
to the maximum and in (di) and (d2) to the minimum of energy vs. fiux curves. 

the resonance moves away from the Fermi energy and 
therefore, the conductance drops as we enter this regime. 
Finally, there are two nonmagnetic regimes, one in which 
the "impurity" level is predominately empty, 64— sp ^ A, 
known as the empty orbital regime, and the correspond- 
ing regime where the dot is doubly occupied. In these 
regimes, the conductance drops toward zero. 

In Fig. 10 the zero-temperature conductance through 
a quantum dot acting as an Anderson impurity is pre- 
sented and compared to exact results of the Bethe ansatz 
approach''^'^^ . To calculate the conductance, Eq. (31) 
was used, with the ground-state energies at ^ = and 
(f) = TV obtained using the variational method presented in 
Appendix B. There arc two variational parameters defin- 
ing the auxiliary Hamiltonian (Bl), one describing the 
effective energy level on the dot and the other renormal- 
izing hoppings into the leads. Two different variational 
basis sets were used in calculations. In the first set. 




Figure 9: The Anderson impurity model realized as a quan- 
tum dot coupled to two leads. The dot is described with the 
energy level £d and the Coulomb energy of a doubly occupied 
level U. ti is the hopping between the dot and leads. 



conductance in the Kondo regime reaches the unitary 
limit of Letting either ed or Sd + U approach the 
Fermi energy so that \ed + U — ep] or \ed — ep] becomes 
comparable with A, wc enter the mixed valence regime 
where the charge fiuctuations on the dot become impor- 
tant. In this regime, the resonance becomes wider and 
merges with the resonance corresponding to Sd or Sd + U 
levels. More important for our discussion is the fact that 
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the basis consisted of wavefunctions (B2). As a result 
of the rotational symmetry in the spin degree of free- 
dom, two of the basis functions may be merged into one. 
Therefore, the basis set consisted of projections of the 
auxiliary Hamiltonian ground state |0) to states with 
empty Pq |6), singly occupied Pi |0) = P^\0) + Pi |6) 
and doubly occupied P2 |0) dot level. In the second ba- 
sis set, wavefunctions Pi^Pq |0> , PqFPi |6), P2VP1 |0) 
and P1VP2 |0) (B7) were added to those of the first set, 
with V = Vl + Vr being the operator describing the hop- 
ping between the dot and the leads (7). For each position 
of the Ed level relative to the Fermi energy, we increased 
the number of sites in the ring until the conductance con- 
verged. The number of sites needed to achieve conver- 
gence (see Fig. 11a) was the lowest in the empty orbital 
regime and the highest (about 1000 for the system shown 
in Fig. 10) in the Kondo regime. This is a consequence 
of Eq. (19) as a narrow resonance related to the Kondo 
resonance appears in the transmission probability of the 
quasiparticle Hamiltonian (47) in the Kondo regime. In 
the mixed valence regime, the width of the resonance be- 
comes comparable to A, which is much larger than the 
Kondo temperature and the convergence is thus faster. 
In the empty orbital regime the resonance moves away 
from the Fermi energy and an even smaller number of 
sites is needed to achieve convergence. 
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Figure 11: (a) Results of conductance calculations using 
Eq. (31) for the system presented in Fig. 10 as the number 
of sites in the ring increases. Note that the convergence is 
the fastest in the empty orbital regime and the slowest in the 
Kondo regime, (b) Finite-size scaling analysis of the same 
results for various values of ea- With black dots, the Bethe 
ansatz values are shown. Energies were calculated using the 
variational method of Appendix B with 7 variational basis 
functions. 



Figure 10: The zero-temperature conductance calculated 
from ground-state energy vs. magnetic flux in a finite ring 
system using the variational method of Appendix B with 3 
and 7 basis functions. For comparison, the exact Bethe ansatz 
result is presented with a dashed line. The system shown in 
Fig. 9 was used, with U = 0.64<;o and = 0.2io. 

Let us return to results shown in Fig. 10. Note that 
extending the variational space from 3 to 7 basis func- 
tions significantly improves the agreement with the exact 
result. The remaining discrepancy at the larger basis set 
can be attributed to the approximate nature of the varia- 
tional method. Another source of error could be the fact 
that the Bethe ansatz solution assumes there is a con- 



stant coupling to an infinitely wide conduction band. In 
our case, the conduction band is formed by the states in a 
tight-binding ring, the coupling to which is not constant. 
However, it is almost constant in the energy interval we 
are interested in, i.e. near the center of the band. In 
order to estimate the effect of the nonconstant coupling 
on the conductance, we calculated the dot occupation 
number within the second order perturbation theory for 
both the case of a constant coupHng and for the case 
of a tight-binding ring. Wc then calculated the conduc- 
tance in each case making use of the Friedel sum rule"*^ . 
The agreement is significantly better than the difference 
between the Bethe ansatz and variational conductance 
curves in Fig. 10. Therefore, we believe that the use of 
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Bethe ansatz results is justified for tfiis particular prob- 
lem. 

In Fig. lib a finite-size scaling analysis of tfie con- 
vergence is presented. Note that for rings with a large 
number of sites N, the error scales approximately as jf. 



C. Double quantum dot 

The next logical step after studying individual quan- 
tum dots is to consider systems of more than one dot. 
Single quantum dots are often regarded as artificial atoms 
because of a similar electronic structure and comparable 
number of electrons in them. By coupling several quan- 
tum dots one is thus creating artificial molecules. Here 
we will not go into detail in describing the physics of such 
systems. Our goal is to compare results of our conduc- 
tance formulae to results of other methods for a double 
quantum dot system presented in Fig. 12. 

In the calculation we again employed the conduc- 
tance formula (31) and calculated the ground-state en- 
ergies with the variational approach of Appendix B with 
the variational basis set (B2). In Fig. 13 the zero- 
temperature conductance for the case where the inter- 
dot and the on-site Coulomb repulsions V and U are of 
the same size, are plotted as a function of the position 
of dot energy levels relative to the Fermi energy for var- 
ious values of the inter-dot hopping matrix element t2. 
The same problem in the particle-hole symmetric case 
was studied recently in Ref.*. The Mat- 
subara Green's function was calculated with the quantum 
Monte Carlo method and the values on discrete frequen- 
cies were extrapolated to obtain the retarded Green's 
function at the Fermi energy. Then Eqs. (40) and (39) 
were used to calculate the zero-temperature conductance. 
The results are presented in Fig. 14, together with the 
conductance calculated within the Hartree-Fock approx- 
imation and results of our method. The agreement with 
the QMC results is excellent, while the Hartree-Fock 
approximation gives a qualitatively wrong conductance 
curve, especially at low values of the inter-dot coupling, 
indicating strong electron-electron correlations in the sys- 
tem. The results of our method for lower values of V are 
also shown in Fig. 14. 




(E^+U/2+V)/t^ 



Figure 13: The zero-temperature conductance of the system 
in Fig. 12 as a function of the position of the dot energy level 
Ed and inter-dot hopping matrix element t2. The remaining 
parameters are U = V = to and ti = 0.5to- 



.... V/t^ = 
__ V/tg = 0.5 

— VltQ = 1,HF 
Vlt^= 1,QMC 




Figure 14: The zero-temperature conductance of the double 
quantum dot system of Fig. 12 at £<£ -1-^-1-^ = as a func- 
tion of the inter-dot hopping matrix element t2 for various 
values of the inter-dot Coulomb interaction V . As a compari- 
son, the Hartree-Fock and quantum Monte Carlo results® are 
presented for V/to = 1. Other parameters are the same as in 
Fig. 13. 




Figure 12: A double quantum dot system. Each of the dots 
with energy level Ed and on-site Coulomb repulsion U is cou- 
pled to a lead with a hopping matrix element ti . The inter-dot 
hopping t2 is also present as is the inter-dot Coulomb repul- 
sion V . 



VI. SUMMARY AND CONCLUSIONS 

We have demonstrated how the zero-temperature con- 
ductance of a sample with electron-electron correlations 
and connected between nonintcracting leads can be de- 
termined. The method is extremely simple and is based 
on several formulae connecting the conductance to per- 
sistent currents in an auxiliary ring system. The conduc- 
tance is determined only from the ground-state energy 
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of an interacting system, while in more traditional ap- 
proaches, one needs to know the Green's function of the 
system. The Green's function approaches are often much 
more general, allowing the treatment of transport at fi- 
nite temperatures and for a finite source-drain voltage ap- 
plied across the sample, which in our method is not possi- 
ble. However, the advantage of the present method comes 
from the fact that the ground-state energy is often rela- 
tively simple to obtain by various numerical approaches, 
including variational methods, and could therefore, for 
zero-temperature problems, be more appropriate. 
Let us summarize the key points of the method: 

(1) The "open" problem of the conductance through a 
sample coupled to semiinfinite leads is mapped on to a 
"closed" problem, namely a ring threaded by a magnetic 
fiux and containing the same correlated electron region. 

(2) For a noninteracting sample, it is shown that the 
zero-temperature conductance can be deduced from the 
variation of the energy of the single-electron level at the 
Fermi energy with the fiux in a large, but finite ring sys- 
tem. The conductance is given with Eq. (18), or with 
three simple formulae Eq. (20), Eq. (21) and Eq. (22). 

(3) Alternatively, the conductance of a noninteract- 
ing system is expressed in terms of the variation of the 
ground-state energy with fiux, Eq. (24). Three additional 
conductance formulae, Eq. (31), Eq. (32) and Eq. (33), 
are derived. 

(4) The method is primarily applicable to correlated 
systems exhibiting Fermi liquid properties at zero tem- 
perature. In order to prove the validity of the method for 
such systems, the concept of Fermi liquid quasiparticles 
is extended to finite, but large systems. The conductance 
formulae give the conductance of a system of noninter- 
acting quasiparticles, which is equal to the conductance 
of the original interacting system. We also proved that 
for such systems, the ground-state energy of a large ring 
system is a universal function of the magnetic fiux, with 
the conductance being the only parameter [Eqs. (29) and 
(30)]. 

(5) The results of our method arc compared to results 
of other approaches for problems such as the transport 
through single and double quantum dots containing in- 
teracting electrons. The comparison shows an excellent 
quantitative agreement with exact Bethe ansatz results 
in the single quantum dot case. The results for a double 
quantum dot system also perfectly match QMC results 
of Ref.8. 

(6) One should additionally point out that in the 
derivation presented in this paper we assumed the in- 
teraction in the leads to be absent. It is clear that this 
assumption is not justified for all systems. The method 
cannot be directly applied to systems where the interac- 
tion in the leads is essential, as are e.g. systems exhibit- 
ing Luttinger Hquid properties. 

(7) The validity of the method is not limited to systems 
that do not break the time-reversal symmetry. A general- 
ization to systems with a broken time-reversal symmetry, 
such as Aharonov-Bohm rings coupled to leads, is possi- 



ble and will be presented elsewhere . 

(8) Another important limitation of the present 
method is the single channel approximation for the leads. 
It might be possible to extend the applicability of the 
method to systems with multi-channel leads by studying 
the infiuence of several magnetic fiuxes that couple dif- 
ferently to separate channels. This way, one might be 
able to probe individual matrix elements of the scatter- 
ing matrix and derive conductance formulae relevant for 
such more complex systems. 

Note added. After the present work was completed 
the authors met R. A. Molina and R. A. Jalabert who 
reported about their recent unpublished work where an 
approach similar to our work is presented^^'^^. 
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Appendix A: FERMI LIQUID IN A FINITE 
SYSTEM 

In Sec. IV we based the proof of the validity of the 
conductance formulae for Fermi liquid systems on the 
assumption that 

E{N,(t>\M + m) = E {N, M) + 

rn 

+ ^e{N,(l>-M-i) + o(N-iy (Al) 

Here E (TV, </>; M m) and E {N, </>; M) are the ground- 
state energies of an interacting N-site ring with fiux 
ct>, containing M + m and M electrons, respectivclj'. 
e (N, (±>: M; i) is a single-electron energy of the ring quasi- 
particle Hamiltonian H {N, M) as defined in Eq. (49), 
with the Fermi energy corresponding to M electrons in 
the system. The index i labels successive single-electron 
energy levels above the Fermi energy. Wq. assume m to 
be finite and N approaching the thermodynamic limit. 
In this Appendix we will give arguments showing that 
the assumption of Eq. (Al) is indeed valid. In Sec. 1 we 
first express the problem in terms of the Green's function 
of the system. In Sec. 2 we study the properties of the 
self-energy due to interaction in finite ring systems and 
then use this results to complete the proof in Sec. 3. 
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Relation to the Green's function 



2. Self-energy due to interaction 



Assume we manage to prove Eq. (Al) for m = 1, i.e. 



E{N,(j);M + 1) = E {N, (p; M) + 
+e(N,^-M-l) + o(N-i'^ . 



(A2) 



Then we can use the same result to relate the energy of a 
system with M+2 electrons to that with M+1 electrons, 

E{N,(j);M + 2) = E {N, (j); M) + 

+e (AT, 0; M; 1) + e (AT, M + 1; 1) + O (iV- i) . (A3) 

Now the matrix elements of quasiparticle Hamiltonians 
H [N. (j)-, M + 1) and H [N. 0: M) differ by an amount of 
the order of To sec this, note that the shift of the 
Fermi energy as an electron is added to the system is of 
the order of j^, producing a shift of the same order in the 
self-energy and it's derivative at the Fermi energy, which 
define the quasiparticle Hamiltonian through Eq. (49). 
As the difference of the Hamiltonians AH is small, we 
can use the first order perturbation theory, 

e{N,(l);M + 1-1) = 

= s{N,(f);M;2) + (N,(j>;M;2 AH Ar,,^;M;2^ = 

(A4) 

In the last step we made use of the fact that the quasi- 
particle Hamiltonians differ only in a finite number of 
sites in and in the vicinity of the central region, and of 
the fact that the amplitude of the quasiparticle single- 
electron wavefunction \N. <i)\ M; 2) is of the order of 

Thus we have proved Eq. (Al) for m = 2 and using the 
same procedure, we can extend the proof to any finite m. 

To complete the proof, we still need to show the va- 
lidity of Eq. (A2). As a first step, consider the Lehmann 
representation of the zero-temperature central region 
Green's function 

G,, (t, t') = -ie {t - t') (0| [dj (t) , 4 it')] |0) (A5) 

of a ring system characterized with N and (j), containing 
M electrons, 

Gji {N,cj)-M;z) = 

^ WM){n\d\\0) ^ 



= e{N, <t>;M-2)+0{N-^) 



z-{E^+^-E^) 

(0|4|n)(n|rfj|0) 
z-{E^-E^f-^) 



(A6) 



The first sum runs over all basis states with M + 1 elec- 
trons, while the second sum runs over the states with 
M — 1 electrons. The difference in the ground-state en- 
ergies of systems with M + 1 and M electrons is evi- 
dently equal to the position of the first (5-peak above the 
Fermi energy in the spectral density corresponding to the 
Green's function. In what follows, we will try to deter- 
mine the energy of this (5-peak. 



To achieve the goal we have set in the previous Section, 
we first need to study the structure of the self-energy due 
to interaction in a finite ring with fiux. Let us again con- 
sider the Lehmann representation (A6) and to be specific, 
limit ourselves to states above the Fermi energy. Intro- 
ducing (/j^ = (0|dj |n) and £„ = E^+'^ - E^' , 
express the Green's function as 



we can 



Gji {N,<i>-M-z) = Y, 



(A7) 



This expression can also be interpreted as a local Green's 
function of a larger noninteracting system, consisting of 
the central region and a bath of noninteracting energy 
levels, the number of which is equal to the number of 
multi-electron states with M + 1 and M — 1 electrons of 
the original interacting system. The self-energy due to 
"hopping out of the central region", which includes both 
the effects of the interaction as well as those due to the 
hopping into the ring, can then be expressed as 



(A8) 



where Vjn are the "hopping matrix elements" between the 
central region and the "bath". Thus wc have shown that, 
as far as the single-electron Green's function is concerned, 
the interacting system can be mapped on a larger, but 
noninteracting system. 

To further clarify the concepts introduced above, we 
calculated the self-energy due to interaction within the 
second-order perturbation theory. Following the calcula- 
tions by Horvatic, Sokcevic and Zlatic^^"^^ for the Ander- 
son model, we sum the second order self-energy diagrams 
shown in Fig. 15, including Hartree and Fock terms into 
the unperturbed Hamiltonian. A lengthy but straight- 
forward calculation, which we do not repeat here, shows 
that one can identify the states n of Eq. (A8) with three 
Hartree- Fock single-electron state indices q = {qi,q2,Q3) 
such that qi and 52 are above the Fermi energy and qs is 
below it (or vice versa), and a spin index s. The "bath" 
energy levels 



(A9) 



and the "hopping matrix elements" related to the self- 
energy for electrons with spin a 



HyUJ^'^^^*, 8 = a, 

(AID) 

are then expressed in terms of the Coulomb interaction 
matrix elements (5) , and the Hartree- Fock single-electron 
energies Sq {N, <p\ M) and the corresponding wavefunc- 
tions \ip'^ {N , (f); M)) . In Fig. 16 the positions of i5-peaks 
in the imaginary part of the self-energy as a function of 
magnetic flux through the ring are plotted. Note that as 
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the flux is varied, the positions of the peaks fluctuate by 
an amount of the order of the single-electron level spac- 
ing which is of the order of jj . The weights of the peaks 
also depend on the flux. A similar behavior is expected 
if higher order processes are also taken into account. 




Figure 15: Second-order self-energy diagrams. 




ynvi 




CO 



Figure 16: Dashed lines show the positions of 5-pcaks (A9) in 
the second-order self-energy corresponding to single-electron 
energy levels of an unperturbed system presented with gray 
lines. 

Finally, let us study the self-cncrgy in the thermody- 
namic limit. We will show that in this case, the self- 
energy is independent of flux and is equal to the self- 
energy of the original, two- lead system, shown in Fig. 1. 
To prove this statement, we consider a self-energy Fcyn- 
man diagram for the central region decoupled from the 
ring, which then is obviously independent of flux. To 
calculate the self-energy for the full system, one should 
insert the self-energy due to hopping into the ring into 
each propagator of the diagram. The self-energy due to 
hopping into the ring is 



e(°)(7V,<A;.) = ^^, (All) 



where Ek are the single-electron energy levels of the ring 
decoupled from the central region and Vki = —ip^tn — 
'^Rtm is the hopping matrix element between site i in the 
central region and the single-electron state k in the ring. 
Vki is expressed in terms of the hopping matrix clement 
tLi between the site i and the ring site L adjacent to 
the central region and the single-electron wavefunction 
"^i = Y^jvqrr sin k at site i, where N is the number of sites 
in the ring. There is also a similar contribution to Vki 



corresponding to the hopping into the right lead. In the 
ring system, the right lead wavefunction can be expressed 
in terms of the left lead one as V'^ = e~*'^V'i with 

k = jf^, if one takes into account the parity of the 
wavefunctions and the effect of the flux. Thus, Eq. (All) 
transforms into 



Sl.f ^ {N; z) 



+ 



2 {tjLtR; 



ie 



+ tjRtLie''^) 



S^.f ^ {N; z) 



N + 1 



^t^^, (A12) 



where {N; z) and ^'ji' {N; z) are the self-energies 
due to hopping into the left and the right leads (each 
with N sites) of the two-lead system. In the third term, 
one can perform the sum over odd n-s and over even n-s 
separately. The sums differ only in sign in the N ^ oo 
limit and therefore, this term vanishes. Therefore, in the 
thermodynamic limit the self-energy due to interaction is 
the same in both two-lead and ring systems. 



3. Proof of Eq. (A2) 

Positions of ^-peaks in the spectral density of the in- 
teracting system correspond to the single-electron energy 
levels of the noninteracting part of the ring Hamiltonian 
coupled to the "bath" according to Eq. (A8). These en- 
ergies can be obtained by solving for zeroes of the deter- 
minant of the inverse of the "local" Green's function 



det 



ujl - H(°) {N, ^ii) - S {N, (I);M;uj + i6) 



= 0. 



(A13) 

What we are going to prove in this Section is that the 
lowest positive solution of this equation corresponds to 
e {N, 0; M ; 1) as required by Eq. (A2). 

We begin by separating the self-energy at frequencies 
close to the Fermi energy into two contributions, one 
(S") due to the "bath" states close to the Fermi en- 
ergy and the other (S') of aU the other states with en- 
ergies which are separated from the expected solution 
of Eq. (A13) by at least an amount of the order of the 
single-electron level spacing A, which is of the order of 
jf. We flrst estimate the second term. Let us divide 
the frequency axis into intervals of width A, each con- 
tributing to the self-energy at |a;| < A an amount given 

by 



e+A 



Pie) 



ds, 



(A14) 



where pji (e) = VjnVniS (e — e„) if the notation of 
Eq. (A8) is used. On average, this contribution corre- 
sponds to that of a system in the thermodynamic limit 
where p (e) is a continuous function and the magnitude 



of each contribution is at most of the order of 



N- 



To 



see this, let us assume p (e) is proportional to (41) 
for all values of £ up to a cutoff of the order of 1. Such 
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an approximation can be considered as the upper limit 
of possible values of p (e) in Fermi liquid systems, if one 
does not take into account the rapidly decreasing tails 
at higher energies, which contribute a negligible amount 
to the self-energy at the Fermi energy. Evaluating the 
above integral, we find that contributions of the intervals 
close to the Fermi energy are of the order of and con- 
tributions of the intervals near the cutoff are of the order 
of ^. Using an analogous procedure, we can also eval- 
uate the derivative of the self-energy close to the Fermi 
energy, with contributions 

rde. (A15) 



P{e) 



In this case, also contributions corresponding to intervals 
close to the Fermi energy are of the order of . If p (e) 
for a finite N is used instead, there are large fluctuations 
about the average value (see the discussion in the pre- 
vious Section) with the amplitude of fluctuations of the 
same order of magnitude as the average value itself. To 
estimate the difference between the finite-system's real 
part of the self-energy (or its derivative) close to the 
Fermi energy and the corresponding quantity for a sys- 
tem in the thermodynamic limit, we note that a sum 
of N quantities, each of them of the order of with a 
standard deviation of the same order of magnitude, has 
a standard deviation of the order of A^" 2 , and therefore, 
we can estimate that for Iwl < A 



S' (iV,0; M-uj + i5) = 
as' {N,(f>;M;uj + iS) 



(A16) 



du! 



i5) 



duj 



O 



(n-'^^ . (A17) 



Note that we do not need to exclude the contribution of 
the interval at the Fermi energy (the one corresponding 
to S") from self-energies in the right hand sides of these 
equations, because the corresponding contributions are 
smaller than 7V^2 as discussed above. Also the errors 
arising from the fact that the right hand sides are evalu- 
ated at a; = instead of at w are only of the order of , 
as discussed in the previous Section. In Fig. 17 a com- 
parison of the self-energies at a finite N and in the ther- 
modynamic limit is presented. Note that in the vicinity 
of the Fermi energy, the real parts of both self-energies 
coincide. 

One can now proceed as in Eqs. (44) and (47), defining 
the rcnormalization matrix Z' (TV. (A: M) = Z + 0{N-^) 
and the quasiparticlc Hamiltonian H' (A'', M) = 
H (iV, (j)] M) + ©(A^^ 2) corresponding to the self-energy 
S'. As shown in the previous Section, the sclf-cncrgics 
of an infinite two-lead system and the corresponding ring 
system are the saine and therefore, the renormalized ma- 
trix elements of H (iV, 0; M) correspond to those of a 
two-lead system. For < A, Eq. (A13) transforms into 




1 "1 

1 (1 ' 1' 

II 




iMll 

^ - J 




_2liiLL 









Figure 17: The (a) real and (b) imaginary parts of the self- 
energy of an interacting system in the thermodynamic limit 



and for N 
Fig. 5. 



400 with 



The system is described in 



det wl-H'(iV,0;M)-S"(iV,0;M;w-t-«(5) =0, 

(A18) 

where the coupling to the remaining "bath" levels has 
been renormahzed as S" = z'i/2S"Z'i/2. Let us for a 
moment neglect this term in Eq. (AI8). As the difference 
A# between Hamiltonians H{N, cf): M) and H' (TV, 4)] M) 
is smah for a large N , one is justified to relate their single- 
electron energy levels using the first order perturbation 
formula 



e'(iV,0;M;l) = 

= e(TV, (t>\M;l)^ (n, (P;M;1 

= e{N,cj);M;l)+o(^N-iy 



AH 



iV,0;M;l 



(A19) 



In the last step we made use of arguments simflar to those 
in deriving Eq. (A4). 

The energy (A19) can acquire an additional shift be- 
cause of the coupling S" . To estimate this shift we first 
note that in the worst case scenario, i.e. when there is a 
single bath energy level which coincides with the quasi- 
particle energy level (A19), the coupling matrix elements 
Vjn in Eq. (A8) must be at most of the order of A^~2 for 
Eq. (41) to be satisfied in the thermodynamic limit. Then 
one can make use of the degenerate first order perturba- 
tion theory, which shows that the quasiparticlc energy 
level is shifted by an additional amount of the order of 
This completes the proof of Eq. (Al). 

As a conclusion, in Fig. 18 we present a comparison of 
the total densities of states for a finite ring interacting 
system within the second order perturbation theory and 
in the quasiparticlc Hamiltonian approximation. Note 
that the states near the Fermi energy are well described 
with the quasiparticlc approximation, whfle the states 
further away from the Fermi energy are split in the in- 
teracting case. Similar results were reported in Ref.'"^. 
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Figure 18: (a) The total density of states of an interacting ring 
system within the second-order perturbation theory, (b) Total 
density of states corresponding to the quasiparticle Hamilto- 
nian. The system is described in Fig. 5. 



Appendix B: VARIATIONAL GROUND-STATE 
ENERGY 

In order to calculate the conductance for interacting 
systems, we first need to devise a robust method that 
would allow us to efficiently calculate the ground-state 
energy of such systems. Note that we need a method that 
would provide us with the energy of a system with a very 
large number of sites in the ring. However, the number 
should still be finite - i.e. we must not perform the calcu- 
lations in the thermodynamic limit. We made use of the 
projection method of Gunnarson and Sclionliammer^''^^^ , 
introduced originally to calculate the ground-state energy 
of the Anderson impurity model, and extended it to treat 
the more general Hamiltonian (1). 

Let us introduce an auxiliary noninteracting Hamilto- 
nian, 

H = Hl + Vl+hP +Vr + Hr, (B1) 

with arbitrary matrix elements describing the hopping 
between the leads and the central region, and the cen- 
tral region itself. Note that these are the same matrix 
elements as the ones being renormalized in the Fermi liq- 
uid quasiparticle Hamiltonian (47). Let us also define a 
Hilbert space spanned by a set of basis functions 



where M is the number of sites in the central region, 
|0) is the ground-state of the auxihary Hamiltonian (Bl) 
containing the same number of electrons as there are in 
the ground state of the original Hamiltonian, and 

Po' = (1 - n^r) (1 - na) , (B3) 
Pi = n,,(l-nis), (B4) 
Pi = Ui^riii (B5) 

are projection operators on unoccupied, singly occupied 
and doubly occupied site i. The original Hamiltonian is 
diagonalized in the reduced basis set introduced above, 

H/soc = ESpa, (B6) 

with H/3a = (V'|9 \H\ il'a) being the matrix elements of 
the Hamiltonian and Spa = {i'i3 \ l^-'a) take into account 
the fact that the basis functions do not form an orthonor- 
mal basis set. The eigcnstate with the lowest energy Ejj 
of this eigenvalue problem is an approximation to the 
ground-state energy of the original Hamiltonian. Varying 
the parameters of the auxiliary Hamiltonian, one can find 
their optimal values which minimize Efj. The solution 
of this minimization problem is the final approximation 
to the ground-state energy. 

Let us consider some simple limits of the problem. In 
the noninteracting case where U = 0, one can choose the 
auxiliary Hamiltonian to be equal to the true Hamilto- 
nian H = H. Then the wavefunction = J2c -P« |0) = 
|0) = |0) is the exact ground-state wavefunction of 
the system. Note that applying the same wavefunction 
ansatz to the interacting case and allowing the matrix 
elements of the auxiliary Hamiltonian to be renormal- 
ized, provides us with the Hartree-Fock solution of the 
problem. Therefore, the variational method introduced 
above always gives the ground-state energy which is lower 
or equal to the corresponding Hartree-Fock ground-state 
energy. In the limit of the central region being decoupled 
from the ring, i.e. Vl = Vn = 0, the variational method 
also yields the exact ground-state energy. To prove this 
statement, let us select the matrix elements of H in such 
a way that in its ground state there are m electrons in 
the central region. Then the basis set (B2) spans the 
full Hilbert space for m electrons in the central region. 
As there is no coupling to the states in the ring, solving 
the eigenvalue problem (B6) provides us with the exact 
ground state of the problem with a constraint of a fixed 
number of electrons in the central region. By varying H, 
all the possible values of m can be tested and the one 
yielding the lowest ground-state energy corresponds to 
the correct ground state of the system. 

The variational wavefunction ansatz can be improved 
by extending the Hilbert space with additional basis func- 
tions, the most promising candidates being of type^^ 

l^lSXjiac) = PfiVxjiaPc |0) , (B7) 

where Vxjia = Vxjic'^^dia + h.c. and A is a lead index, 
i.e. either L or R. On the other hand, as the size of the 
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Hilbert space increases exponentially with the number 
of sites in the central region, it might be convenient to 
limit the basis set to the states obtained by projecting 
to the central region's many body states between which 
fluctuations are possible. 

Finally, we state some technical details concerning the 
evaluation of -ff/Sa and S^a- It is convenient to express 
these matrix elements only in terms of quantities related 
to the central region and the neighboring sites in the 
leads. As 



Spec = (6 \PpPoc\ 0) = (0 |P„| 0) 



(B8) 



the scalar products between the basis functions arc evi- 
dently expressed with the central region quantities. The 
matrix elements of the Hamiltonian can be expressed as 



H^o. = {0\PpHP^\Q) 
'q HPpPo. 

: HP^Pc 
ESfSc 



- 



0) + (0|P^i/P„|0>- 

+ {0\Pl3{VL+Hc + VR)Pa,\0)- 

(Vl + H^c^ + Vr) P« I O) 60ac, (B9) 



where E is the ground-state energy of the auxihary 
Hamiltonian H. In the second and the third term we 
made use of the fact that lead Hamiltonians Hi, and Hr 
commute with the central region projectors and there- 
fore, they cancel out. Again, we succeeded in expressing 
the matrix elements in terms of central region quantities 
together with quantities related to the neighboring sites 



in the leads. Similar results are obtained if the extended 
basis set of Eq. (B7) is used. The matrix elements in 
Eqs. (B8) and (B9) need to be calculated in a noninter- 
acting state. Therefore, we can make use of the Wick's 
theorem to decompose the expressions into two-operator 

averages of type ^6 dj^di^ 0^. As a huge number of 

terms is generated in this procedure, the decomposition 
was performed automatically by symbolic manipulation 
of operators. The ground-state energy of the auxiliary 
Hamiltonian and the two-operator averages can be ex- 
pressed in terms of the single-electron energies £fe and 
wavefunctions 



\(p'') of H as 



E 



2 Yl ^k, 



k occ. 



(BIO) 
(Bll) 



k occ. 



The sums run only over the single-electron states occu- 
pied in the ground state |0). The eigenvalues ik were 
calculated in a basis in which the Hamiltonian matrix 
is banded, i.e. linear combinations of local basis func- 
tions corresponding to the left lead and right lead sites 
were introduced to "move" the hopping matrix elements 
in corners of the matrix close to the diagonal. For each 
eigenvalue, only the components of the eigenvector re- 
lated to the central region and neighboring sites were 
calculated, again taking the special structure of the ma- 
trix into account. The procedure used scales with the 
number of sites in the ring as O (iV^), which ahows one 
to treat systems with up to 10000 sites in the ring. 
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